UNIBAS-MATH 10/96 



The WKB Approximation Without Divergences 



D. Cocolicchio^ 1,2 ) and M. Viggiano^ 1 ) 

1 ) Dipartimento di Matematica, Univ. Basilicata, Potenza, Italy 

Via N. Sauro, 85, 85100 Potenza, Italy 

2 ) Istituto Nazionale di Fisica Nucleare, Sezione di Milano, Italy 

Via G. Celoria 16, 20133 Milano, Italy 



ABSTRACT 



In this paper, the WKB approximation to the scattering problem is developed with- 
out the divergences which usually appear at the classical turning points. A detailed 
procedure of complexification is shown to generate results identical to the usual WKB 
prescription but without the cumbersome connection formulas. 



I. Introduction 



In general an exact solution for many quantum mechanical problems is unavoidable and 
one is forced to resort to some type of approximate technique. One of the most useful of 
those methods is the semiclassical or Wentzel-Kramers-Brillouin (WKB) approach [ 1 ]. 
The major shortcoming of the semiclassical WKB approximation of solving the wave 
equation is its divergence at the classical turning points. Presently, available regular- 
ization schemes are accurate, but rather complicated. Although these methods sharpen 
the threshold effects, nevertheless exact solutions to the stationary wave equation 



d z ip 
dx 2 



+ k 2 (x)ij = (1) 



(with the local wavenumber k 2 (x) = 2m[E — V(x)]/h 2 ) cannot be found in most prob- 
lems which involve one-dimensional potential V(x) [2]. In the WKB approach, the 
wave function ifj(x) is supposed to be represented by 

tP{x) = Ae ia{x)/h (2) 

which converts the linear, time-independent Schrodinger equation for ip(x) into the 
non-linear Riccati equation for the function cr(x) 



The WKB approximation consists in expanding a as a power series in h: 



% ( T% \ ^ 

a(x) ~ a (x) + -ai (x) + f - J a 2 (x) + o(h 3 ) (4) 



and substituting this expression in the relative differential equation, whose coefficients 
can be generated recursively [ 3 ] . The following conditions can then be obtained: 

fdaoMV = 2m[E-V(x)]=p 2 (x) 

(5) 

/ dcr (x) \ ( dai(x) \ , if d 2 a (x) \ _ ^ 
\ dx J\dx J "r 2 ^ dx 2 J U 
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whose solutions are given by 



<Jo(x) = ± / p(x') dx' 

Jxo (6) 
^ {x)= -2^y = -2HxJ^ -i(-) = -In 



where the prime denotes differentiation with respect to x and xq is an arbitrary point. 
? le, 

V(xc ) = -E, which makes the first-order WKB solution 



The leading connection term <j\{x) diverges at the classical turning points x^ where 



ip(x) ~ ^ = |c + exp ^ y p(x') cte'^ + C_ exp (^~^ J P( x ') dx ^j | ( 7 ) 



divergent in these points. In the classical limit, this divergence is understandable, since 
a classical particle has zero velocity at these turning points. This divergence at the 
turning points is a severe limitation on the usefulness of the WKB approximation for 
quantum mechanics, since there is no divergence in the exact wave function. However 
away from the turning points, the WKB solution gives a good description of wave 
functions, especially in the semiclassical limit of large quantum numbers n. Far from 
the turning points, the behaviour of the semiclassical solution 




changes drastically in accord with the wave-number 

1 1 

k(x) = -p(x) = -^/2m[E -V{x)} . (9) 

It It 

An oscillatory behaviour is produced by the solution corresponding to the classically 
allowed region E > V max (x), where k(x) is real: 



V>(:r) = -^L= sin ( f k(x') dx' + - X (10) 
y/p(x) [Jxo 4 J 

whereas in the classically forbidden region E < V m i n (x) where the wave-number k(x) = 
i(3(x) with 

P(x) = \y/2m[V{x)-E] > , (11) 

/ 1 



becomes purely imaginary, the general semiclassical solution is exponentially decrescent: 



1>{x) = -/== |c + exp (- J* P(x') dx'^j + C_ exp QT (3(x') dx'^j J . (12) 

This semiclassical behaviour is valid only asymptotically. The regions near the turning 
points should be treated separately. This leads to the fragmentation of the x axis into 
several regions with connection formulas for going through the turning points. Such 
patchwork for the semiclassical wave functions leads to the familiar lowest-order WKB 
energy quantization condition 

J \/lm\E n - V(x)] dx=\n + -J nh (13) 

or to more complicated conditions if higher orders corrections in h are kept [ 4 ] . On 
the other side, continous connection formulas giving finite wave functions also at the 
turning points have also been developed in the realm of uniform approximation [ 5 ] . 
In the next section we recover the connection formulas of the usual WKB procedure, 
then in a later section we establish that this prescription can be generalized with a 
reformulation based on the analytic continuation into the complex momentum variable 
and contour integrations. 



II. Semiclassical Approach to Barrier Penetration 



In the case of scattering problems, there are two independent solutions which are usually 
called incoming and outgoing waves, respectively. These waves become free-particle 
waves in the asymptotic region. Classically a particle (incident from the left with 
energy E < V max {x)) is completely reflected from the potential region at the (left- 
hand) classically turning point a, defined by 

V(a) = E . (14) 

However, quantum mechanically the particle can "tunnel" through the barrier and find 
itself on the right-hand side. In the case of ordinary barrier penetration V m i n < E < 
Vmax there exist two classical turning points on the real axis xi^ = a < = b. The 
WKB approximation is valid where transmission dominates over reflection and, in the 
semiclassical limit, the probability of tunneling is given by T = e~ 2cr * where 

a*= (3(x)dx=- y/2m[V(x) - E] dx (15) 

J a M J a 
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and the probability of being reflected is correspondingly reduced from its classical value 
of unity to R = 1 — T. The problem of the breakdown of the WKB solution near 
the turning points can be overcome with the usual procedure to consider the linear 
approximation of the potential 

V(x) = V(a)+fi(x-a) (16) 

with 
and 

p(x) = >/2m(x(x — a) . (18) 

Thus, the connection formulas, which relate oscillatory and exponential behaviour of the 
wavefunction forms on the opposite sides of a classical turning point, can be matched 
by considering the differential equation 

l! " z\ ili{z) = (19) 



dz 2 

1 /3 

once we set z = (^m/i/h 2 ) (x — a). The complex method for treating classical turning 
points provides a powerful description to the "connection problem" by using the exact 
solution of the Schrodinger's equation 

ip(x) = |c+(a;)exp ^w(x ,x)j +C_(x)exp ^-^w(x , x)j j (20) 

where 



w{xq,x)= I p(x')dx'. (21) 

J x 

The usual WKB method, by its very nature, cannot take into systematic account the 
modifying effects of the multiply reflections between the classical turning points. Even 
if the values of the WKB multipliers C±(x) are known on a given wide region, the 
complex method, while providing useful information about the solutions, it results in- 
efficient to give all the details of the wave function in the neighbourhood of the turning 
points. A powerful technique to overcome this deficiency consists in deriving the asymp- 
totic solutions of the differential Airy equation (19) calculated in the stationary phase 
approximation. Its solution can be represented by the Laplace integral 

i/>(x) = A J e zt ~ t3 / 3 dt (22) 
c 
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where the curve C is taken so that the integrand vanishes at the limit of integration. 
Since the integrand vanishes exponentially in the following interval 



\argt\ < | (23) 

| < argt < (24) 
7 3 

-7r <argt< -n (25) 

the curves C\, C2, C3 are all allowed and yield three different solutions and two each are 
independent. The integration along imaginary axis C\ yields the so called Airy function 

Ailz) = J-/e*p Lt - t\ dt = cos Lt + *\ M (26) 

Cl 

whereas the integration along a complementary curve C2 yields a further independent 
solution 



Bi(z) = ^7 y exp ^ - ^ eft = y sin ^zt + ^ d£ (27) 

c 2 



whose well-known asymptotic behaviours are given by 

^>~^xp(4~ 3/2 ) <28) 

Similarly in the case of Bi(z) we have 

si <-~> ~ ^7I ex P (r S/2 ) ( 3 °) 

Bi( - 2 )._L C0S (^ + ^ . (31) 

These two Airy functions are independent solutions of the Eq.(19). If we introduce the 
variable 

r=\{-zf /2 (32) 
6 



and make the transformation 



= (-^) 1/2 0(r) (33) 

Eq.(19) takes the Bessel form 

^^ + fr'-JW) = (34) 



dr 2 dr V 9 , 
which let us represent Airy functions in terms of the following Bessel functions 

^)=x{^(H-^(H} (35) 
^>=^{MH +/i/s (H} ■ < 3 «> 

The formal connection formulas can be established by the analytic features of the po- 
tential barriers near a classical turning point to match approximate solutions across the 
boundaries. The derivation and application of the connection formulas are both non 
trivial and fraught with pitfalls associated with the existence of exponentially large and 
exponentially small components of the wave function, in the classically forbidden region. 
However in spite of these difficulties, in the case z > (dV/dx > 0) we have the local 
wavenumber in the forbidden region 

P{x) = ly/2m[V(x)-E] = IV^M 

x - a) (37) 

ft It 



so that 

(3(x')dx' = ^z 3 / 2 . (38) 



f 



Similarly for z < we have in the allowed region 

k(x) = \^2m[E - V(x)] = \ ^ -2m[i{x - a) (39) 

It It 

and 

k{x')dx' = -{-zf' 2 . (40) 
x 3 

Finally, since sin (0 + f ) = cos (</> — j) we can derived the following matching expres- 
sion 



=cos 



Since cos (0 + f ) = — sin (</> — j) we find 



sin | / k(x') dx' — - ) < — > exp ( [ d(x') dx' 

Hx) \J X K ' 4 ; y \j a I K ' 



(42) 



^/W> \J X 4; v /^y 

The connection relations for the case of a decreasing potential (dV /dx < 0) are given 
analogously 



ir xp (- r ^ di< ) "vfe™ 8 (r k{x,) dx '-^) 



and 



(43) 



b exp (f ^ "t^T (f fc(a/) ^ " 1) • (44) 



v^M VA y v^( 

With these connection formulae, it is straightforward to determine the transmission 
coefficient T with the knowledge of the semiclassical wave function 

1>(z) = -/== (C+e*- + CLe"*-) (45) 



and its asymptotic limit 



</>fx) ~ <! ^ Iin{ - X) fora; « a Me) 
' ipniout(x) for o< x v ; 



in the region to the far left and right of the barrier, which are given by 

^iin(x) = - --^==[exp ^- J (3(x)dx^j sin {^J k(x')dx - + 
+4iexp(y (3{x)dx j sin (^j k{x')dx' + ^ 



(48) 



and 

ipniout(x) = 7== exp I z / k(x')dx' - i- 
Jkix) V A 



xAM V A ' 4, 

let us make an immediate determination of the transmission coefficient T. This prob- 
ability of tunneling is defined by means of the ratio between the probability current 
density 

j = Re[ —if) — = -1— V ~r - V'-j— (49) 
\zm ax / 2m \ ax ax / 
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of the transmitted and the incident waves: 



T = 



Jin 
ji 



(50) 



where 



3 m = 4\B 



2 n 



3! = A\B\ 2 (e a * + \e~ a *) 2 - (e a * 



A 

rn 



(51) 



with 



a*= 0(x) dx = — I y / 2m[V(x) - E] dx . 

J a Tl J a 



(52) 



Thus the WKB techniques show that in the semiclassical limit, the probability of tun- 
neling is given by 



Jin 



-2(7 * 



"2(7 * 



3i {l + \e-^)' 



(53) 



that is valid under the assumption that exp(— 2a*) 1. Therefore the probability of 
being reflected is correspondingly reduced to R = 1 — T. 



III. Effective Semiclassical Approximation to Barrier Penetration 

In the complex method we previously analyzed, no attempt is made to clarify the 
contributions of the multipliers C±(x) which indeed are ^-dependent and associated 
to the modifying effects of the possibility of internal reflections inside the barrier, in 
the sense that the particle can travel from turning point a to turning b in a several n 
number of ways. A systematic account of multiply reflected contributions in the limit 
of a continuous potential is associated with the following coupled first-order equations 

C'±(x) = a' 1 (x)C^(x) exp ^^w(x , x)^j (54) 

where the prime refers to the derivation with respect to x [6]. These coupled equa- 
tions are formally equivalent to the Schrodinger's Eq.(l). The different solutions of 
the Schrodinger's equation are obtained by applying different conditions to C±(x). In 
particular, imposing that there is no reflected wave far beyond the potential barrier and 
assuming that the incident wave has unit intensity, we obtain the scattering solution, 
after a trivial change of normalization, 

C+(x) = 1 + X! dx'r(x')C-(x')ex.p [-fw(x ,x')] 

(55) 

C-(x) = -C°°dx'r(x')C + (x')e W [+f w ( Xo , x ')} . 
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Wherever the differential reflection coefficient 



r(x) = —a[(x) = 



p'(x) 
2p(x) 



(56) 



may be set equal to zero, a constant value for the C±(x) are obtained, as we mentioned 
above. Therefore, the precise value of the reflection coefficient may be found with the 
help of successive integration by parts 



/+oo 
dx r(x)C+(x) exp 
-oo 
+00 

dx r(x) 



2i . 
+—w(x ,x) 



dx' r(x')C- (V) exp 



+00 



dxr(x) exp 



2i . 
+—w(x ,x) 



~w(xq,x') 



exp 



^■w(xo,x) 



(57) 



In this case, however, the main contributions to the integral cannot be selected easily 
because they depend on the analytic properties of the function p(x) and, in the last 
view, on the type of the singularities of the potential. In fact, the Schrodinger Eq.(l) 
takes the form 



d 2 (j)(w) 
dw 2 



+ 



1 



+ V(w) 



(f>(w) = 



(58) 



where we change the variable to the phase w of the exponent in Eq.(54) and make the 
transformation 



<Kx) 



p(x) 



ip(x) 



being V determined by p through the relation 



V 



3(p'f-2pp" a^' + cn' 2 



4 p 4 



(TO' 



(59) 



p(x) 



02 



(60) 



Here prime denotes differentiation with respect to x. Nevertheless, we may stress that 
the most important quantity results the phase w of the exponent. The subtleties in- 
volved in the evaluation of the precise value of the reflection coefficient stem on the 
convergent expansion of the total wave amplitude which can be given in momentum 
space as 



iP(k 0l k)=<f>(k (h k)+ / dk , Go(ko,k')V(k ,k , )<j ) (ko,k , ) + 



+ / dk' 



dk' G (k , k)V(k , k')G (k', k")V(k', k")<f>(k', k") 



(61) 



+ ... 
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where <f> represents the free wave solution and 



G (k) 



2tt [\-{hk) 2 



(62) 



is its relative propagator. Both <j) and Go are obtained assuming that the perturbation 
V vanishes in Eq.(58). Using these asymptotic solutions, we derive that the exact 
solution of reflection coefficient is perturbatively given by 



Such results would need to take account of all the singularities of the perturbation V 
and indeed of p(x). Thus any value of R will not in general be single valued because of 
the branch points at the real turning points x£ where p(x^) = 0, unless we adopt the 
convention of dividing the complex Gauss plane into two Riemann sheets as shown in 



Concluding Remarks 

The propagator (Green's function) technique in the solution of problems in non rel- 
ativistic Quantum Mechanics becomes relevant in calculating explicitly the reflection 
coefficient in the case of barrier scattering represented by an analytic function V(x). In 
the case of ordinary barrier penetration E < V max , k(x), defined in Eq.(9), is in general 
two-sheeted, with two branch turning points Xc , X c located on the real axis and we 
may choose the defining branch cut to connect them (Fig. 1). An alternative approach 
to the semiclassical approximation allows a very appealing picture of the transmission 
coefficient and generates results identical to the usual WKB prescription, but without 
the cumbersome connection formulas. This method consists in the analytic continuation 
into the complex momentum variable and contour integrations, wherein it is permitted. 
Such a method of complex variable in modern theoretical physics is extensively adopted 
to clarify the concepts of analyticity in S-matrix theory. It results particularly suitable 




(63) 




(64) 



Fig.(l). 
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for discussing the problem of one dimensional barrier penetration if we reconsider that 
the generic propagation from x\ to X2 far to the left and right of the barrier respectively 
in opposite sides can be considered as occurring in three successive steps according to 
the decomposition of the outstanding integral 




The complexification of the problem generates a more intuitive picture of the physics of 
the process and offers an alternative technique to avoid singularities which occur in the 
standard WKB. The propagator can be, then, expressed in terms of contour integrals 
which connects the initial and final points Xi n , Xfi n , both of which are located on the 
real axis far to the left of the barrier. Furthermore we assert that if the incident wave is 
initially located in Xi n far to the left of the barrier, then the reflected wave is given by the 
analytic continuation of the functions involved, evaluated on the other side of the cut. 
It is important to note that they are on different sheets so that any independent contour 
of integration has to pass through the cut. Of course, the singularities of the function 
V in the complex plane may involve other branch points, but we may discard them 
here. Anyway, the singularities of V are related to the singularities of the function 
p(x). Clearly, the multiple integrals resulting in Eq.(63) correspond to the effect of 
multiple reflections. The contributions of the only once-reflected waves to the reflection 
coefficient are given by 

/+oo 
dxr(x)e 2iw{x °' x)/h . (66) 
-oo 

Actually, one could include also contributions from contours which loop the branch cut 
several times. These higher order terms are, in general, unreliable although they are 
expected exponentially much smaller than the polynomial corrections which lie outside 
the usual semiclassical WKB approximation. 
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